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The spectra of Pu chalcogenides and pnictides are computed with LDA+DMFT and interpreted 
with the aid of valence histograms and slave-boson calculations. We find the chalcogenides are 
mixed- valent (n/ = 5.2) materials with a strongly T-dependent low-energy density of states and a 
triplet of quasiparticle peaks below the Fermi level. Furthermore, we predict a doublet of reflected 
peaks above the Fermi level. In the pnictides, the raising of /® states relative to suppresses 
valence fluctuations, resulting in integral- valent (n/ = 5.0) local moment metals. 

PACS numbers: 71.27.+a, 75.30.Mb 



The stark contrast in behavior between the plutonium 
monochalcogenides and monopnictides is a longstand- 
ing issue in strongly-correlated physics. The pnictides 
(PuSb, PuAs, PuP) are comparatively simple metals [1] 
with embedded /-moments arising from trivalent Pu ions 
which order in the range Ocw — 85 to 126 K [2]. In con- 
trast, the chalcogenides (PuTe, PuSe, PuS) exhibit seem- 
ingly contradictory behavior: they have a large room 
temperature specific heat [3] , yet the resistivity indicates 
a small gap [4, 5, 6]. Also, their lattice constant rules 
out the full-shell divalent Pu state, yet the susceptibil- 
ity shows no evidence of Curie- Weiss behavior [2]. Fur- 
thermore, photoemission observes a triplet of peaks (the 
"photoemission triplet") near the Fermi level [7, 8, 9, 10] 
whose origin is still hotly debated [7, 8]. 

The contrast between the Pu chalcogenides and pnic- 
tides exemplifies the view that the Pu 5/ electrons sit at 
the edge of a localization-delocalization transition, where 
small changes in their electronic environment can drive 
a transition to itinerancy or localization, thus posing a 
major challenge to electronic structure methods. Theo- 
retical studies of the chalcogenides within LDA [11, 12] 
predict a metal and do not account for the photoe- 
mission triplet. Methods treating correlations beyond 
LDA have improved the situation, but cannot fully in- 
tegrate the available experimental data within a sin- 
gle theory. Non-charge-self-consistent LDA-I-DMFT with 
FLEX predicts metallic behavior and misses the photoe- 
mission triplet [13]. Adding charge-self-consistency [14] 
opens a gap, but still misses the photoemission triplet. 
LDA-hDMFT with either ED in a small Hilbert space [15] 
or Hubbard-I [16] as the impurity solver describes the 
photoemission triplet but fails to explain the resistivity. 

In this Letter, we elucidate the mechanism driving the 
electronic trends between the pnictides and chalcogenides 
within a single framework. We find that the chalco- 
genides are mixed-valent compounds where valence fluc- 
tuations combine with the underlying Pu atomic mul- 
tiplet structure to drive the formation of a multiplet 
of many-body quasiparticle peaks ( "quasiparticle mul- 
tiplets") which correspond to the observed photoemis- 
sion triplet. These heavy quasiparticles strongly affect 



the density of states at the Fermi level as a function of 
temperature, corroborating the gap-like resistivity and 
large specific heat at room temperature. Using analytic 
methods, we provide a description of the quasiparticle 
multiplet formation and their coexistence with the de- 
velopment of a gap. In contrast, the chemistry of the 
pnictides shifts the atomic multiplet energies, rendering 
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FIG. 1: Computed spectrum for (a) PuTe and (b) PuSb com- 
pared with photoemission [10] . A triplet of peaks are present 
in PuTe, as well as a predicted reflected doublet of peaks (ar- 
rows), while neither appear in PuSb. Application of broad- 
ening (40 meV) has blurred the gap in PuTe (see Fig. 2). 
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FIG. 2; Detail of PuTe spectrum near Fermi level, showing 
development of gap and formation of main quasiparticle peak 
with decreasing temperature. 



valence fluctuations too costly, thereby localizing the / 
electrons. 

We use LDA+DMFT [17, 18] with OCA as the impu- 
rity solver to model the chalcogenides and pnictides, tak- 
ing PuTe and PuSb as representatives of the two groups 
due to the special attention [1, 4, 6, 7, 10] accorded to 
them in the available experimental data. In our calcula- 
tions, we use the projective orthogonalized LMTO basis 
set [19]. We use U = 4.5 eV for the Coulomb interaction, 
consistent with previous work [20, 21, 22]. The Slater 
integrals F'^ = 6.1 eV, F'^ = 4.1 eV and F^ = 3.0 eV are 
calculated using Cowan's atomic structure code [23] and 
reduced by 30 % to account for screening. 

In Fig. 1, we show the computed spectral functions for 
PuTe and PuSb, resolution-broadened by 40 meV and 
overlaid with experimental photoemission data [10]. The 
calculations clearly corroborate the presence of the pho- 
toemission triplet in PuTe at the correct energies, and 
their absence in PuSb. Furthermore, we predict the exis- 
tence of a doublet of peaks (arrows) in PuTe at reflected 
energies about the Fermi level. The strong temperature 
dependence of all five peaks indicates they are quasipar- 
ticle resonances. Examining the main quasiparticle peak 
at the Fermi level (Fig. 2), we find it is composed of 
heavily renormalized quasiparticles with Z w 0.1, giving 
a greatly enhanced specific heat. Additionally, the peak 
sharpens with decreasing temperature, considerably re- 
ducing the density of states at the Fermi level, leading to 
the formation of a gap and the observed temperature de- 
pendence in the specific heat [3]. Together, the reduction 
of Fermi level density and heavy renormalization explain 
how a gap-like resistivity can coexist with a large specific 
heat coefficient at room temperatures. 

A useful way to analyze the Pu atomic environment 
is to quantify the amount of time the /-electrons spend 
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FIG. 3: Valence histograms obtained by projection of DMFT 
solution (r = 60 K) onto Pu atomic eigenstates, plotted with 
energies relative to lowest-energy atomic state. The height of 
each bar represents the percentage of time the atom spends 
in each configuration. PuTe is strongly mixed-valent due to 
the small energy cost (~ 0.5 eV) of valence fluctuations from 
f'' to while PuSb is integral valent due to the large cost of 
fluctuations to both and In PuTe, additional fluctua- 
tions to atomic multiplets ^ H-j/2 and ®-f5/2 are crucial to the 
creation of the photoemission triplet. All quantum numbers 
are gathered into a single index m or n used in the slave-boson 
calculation. Approximate term symbols are given although L 
and S strictly are not good quantum numbers. 



in each atomic conflguration as they fluctuate between 
the atom and conduction band. To this end, we project 
the DMFT ground state \Q) onto the Pu /-electron 
atomic eigenstates, resulting in the probabilities P„i — 
Z^'^ {Vt\X,nm\^) , where Z = I]m(^l^™ml^) is the nor- 
malization and Xmrn is the Hubbard operator which 
projects onto the mth atomic eigenstate [22]. Plotting 
P,„ against the atomic energies gives a valence histogram 
(Fig. 3) which graphically represent the relative weights 
of the atomic configurations comprising |f2). The / va- 
lence can then be defined by (uf) = J2m-^rnnm, where 
Urn, is the number of electrons in the mth state. 

The histograms show that the Pu atom is restricted to 
just one or two valences in both compounds. In fact, the 
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atom mostly exists in a single configuration 
tall blue bar), which we loosely call the "ground state". 
However, PuTe differs from PuSb in that its electrons 
have an over 20 % probability to fluctuate to the low- 
est state CFq, leftmost red bar) due to this state's 
small 0.5 eV separation from the "ground state" . The re- 
sulting mixed- valent {{nf) = 5.2) Pu atom strongly sug- 
gests a Kondo-like, and thus nonmagnetic, ground state 
in PuTe. Additionally, the proximity in energy of the 
next two higher states at 0.5 eV and 

at 0.9 eV) renders these multiplets accessible to valence 
fluctuations, which will play a role in generating the pho- 
toemission triplet. In contrast, the Pu atom is integral- 
valent ((n/) = 5.0) in PuSb. Within LDA, the j = 5/2 
/-bands are higher in energy in PuSb, raising the energy 
of the states relative to the The resulting 1.5 eV 
gap locks Pu into the lowest state and PuSb remains 
a local moment metal. 

To gain additional insight into the LDA-I-DMFT solu- 
tion, we construct a Hamiltonian for the DMFT quantum 
impurity. The histograms indicate we only need to keep 
two valences in the atomic Hilbert space for a low-energy 
model, 



(1) 



where the auxiliary fermions /m|0) — |m; f^) and bosons 
&JjO) = 1"-;/^) create the atomic eigenstates, and E^^ 
and E^-^ are the corresponding atomic eigenenergies. The 
Hamiltonian is supplemented by the constraint Q — 
/m/m + E„ = 1 , in the Same spirit as the slave- 
boson construction [24, 25]. The atom hybridizes with an 
auxiliary conduction bath, 



ka 



+ V(14a4cfea+ll.C.), (2) 



ka 



where c?Jj creates an electron in the ath atomic crystal 
field basis and k is the dispersion of the conduction bath. 
Since we work in the atomic eigenbasis, we eliminate d) 
in favor of the auxiliary particles by expanding d^, = 
bi{F°'^)nmfni, whcrc (F"'f)„„ = {niSjm) are the matrix 
elements of the physical electron creation operator. 

This model is equivalent to the slave-boson treatment 
of the multi-orbital Anderson impurity model [26, 27], 
so we can compute the mean-field solution and fluctu- 
ations. At the mean-field level, we replace the bosonic 
operators by their averages, (bl)^ = {bn)^ = z„ which 
are the probabilities of the atomic states (red bars in 
Fig. 3). Then, the physical Green's function is 



where the auxiliary / propagator and hybridization are 

G^iuj):^]^^-^ = {iu ~ £^,4 ~ ^)'5m'm + i^m'rn SgnU, (4) 

^m'ni — ^ ^ ^m'n'i-^ ^)nmy/ ^n' ^n^aa ; (5) 
n'na 

and the hybridization Aaa' is approximated as an energy- 
independent constant. Here, A is the Lagrange multi- 
plier used to maintain (Q) ~ 1. The crucial minus sign 
GI^i ^{—iiS) arises because the propagation of a physical 
electron ^ {d{T)S) corresponds to an /-hole. 

In the Kondo regime, the mean-field equations give 
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where z = Zn is the total /^ probability and {E^) = 
z^^ J2n ^n^n is the weighted average of the /^ energy 
levels. We ignored the off-diagonal components of gI^,^, 
which are negligible compared to the diagonal compo- 
nents when \E^, — EI^\ 3> A^'m- In the Kondo regime, 
A ~ ~F!q , pinning the lowest / propagator near the 
Fermi level. 

For PuTe, we explicitly evaluate the sum in (3) to de- 
termine the origin of the photoemission triplet. Keeping 
just zq w 0.20 since the remaining z„ are negligible (see 
Fig. 3), we find that only three matrix elements have 



significant weight: Fqq"^ , FqI^ and Fji^- Again ignor- 
ing off-diagonal terms, we find for the physical Green's 
function 

G5/2(^) = zo\F^I'\'gU-u^) + zo|Fo2?G^2(-^), (7) 
Gy^ico) - zo|i^oY?Gfi(-c.). (8) 

The selection rules contained in the matrix elements de- 
termine the spin-orbit structure of the spectrum: G5/2 
contains the main peak at the Fermi level and the weaker 
peak at —0.9 eV (labeled by brackets in Fig. 1), while 
Gj/2 contributes the third peak at —0.5 eV. Plugging 
in £) = 1.0 eV and A5/2 = A7/2 = 0.04 eV into (6) 
gives Tk ~ 500 K in PuTe, so the peaks have sufficient 
width to be seen in photoemission. Numerical solution of 
the mean-field equations with a small hybridization gap, 
Aaa{i^) — Aaa[&{oJ — Eg) + 9{—uj — Eg)], coufirms that 
for 2Eg < Tk, the resonances are not destroyed. Thus, in 
PuTe, valence fluctuations create three Kondo peaks in 
the DMFT quantum impurity with spacings determined 
by the underlying atomic multiplets, corresponding to 
the quasiparticle triplet in the lattice. 

Proceeding to PuSb, the decrement in valence allows 
Pu to fully transfer three electrons to the pnictogen and 
exist purely in a trivalent state. Energetically, this is 
accomplished by raising the /^ multiplets with respect to 
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from 0.5 eV to 2.0 eV. This exponentially suppresses the 
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ground state excited singlet excited singlet pliotoelectron 
singlet I GS) at 0.5 eV at 0.9 eV state dt|GS) 

FIG. 4: Valence histograms of states responsible for the 
predicted doublet of peaks above the Fermi level in PuTe. 
Histograms (a)-(c) depict the ground state and two excited 
Kondo singlets respectively. All three state have significant 
/® weight (arrows) due to strong valence fluctuations. His- 
togram (d) is the photoelectron state d^|GS) generated by the 
addition of an electron to the ground state. Since destroys 
fm and creates 6^, the result is a state which is entirely f'' 
in character. The large admixture in the excited singlets 
means they will overlap strongly with the photoelectron state, 
giving rise to the doublet of peaks at 0.5 eV and 0.9 eV. 



Kondo temperature of PuSb to under 1 K, eliminating 
the Kondo peaks, localizing the /-electrons and allowing 
magnetically-ordered states at low temperature. 

To explain the Hubbard bands and reflected doublet of 
peaks above the Fermi level, we compute corrections to 
the mean- field solution. These corrections show that the 
bare atomic multiplets generate Hubbard bands, which 
are too broad to be seen in Fig. 1. We emphasize that the 
photoemission triplet is not directly attributed to atomic 
multiplets, but rather to quasiparticles. Additionally, the 
corrections to mean-field show that the reflected doublet 
of peaks arises from overlap of the ground state Kondo 
singlet with two excited Kondo singlets (Fig. 4). While 
the ground state is primarily a singlet formed between 
/o and a conduction electron, the two excited states are 
singlets formed with fi and /2 in place of /q. Since the 
atomic multiplets /i and /2 lie at energies 0.5 eV and 
0.9 eV above /o, the two excited Kondo singlets lie at 
these energies as well. To generate the spectrum, an 
electron is first added to the ground state singlet to create 
a photoelectron state (Fig. 4d), which is overlapped with 
the excited singlets. The overlap is non-zero since all 
states have an component due to valence fluctuations, 
resulting in a doublet of reflected peaks. 

To summarize, our LDA-I-DMFT calculations account 
for the complex trends observed in experiment across 
the Pu chalcogenides and pnictides. These include the 
formation of a low-temperature gap at the Fermi level 
accompanied by quasiparticle multiplets in the photoe- 
mission spectrum. Our theory elucidates the mechanism 
for the emergence of the quasiparticle multiplet excita- 
tions. These excitations represent the remnants of the 
atomic structure in the low-energy spectra, which is en- 
tirely described in terms of quasiparticles. The spectra 
and /-occupancy confirm that the chalcogenides are cor- 
related low-carrier materials in the mixed-valent regime. 
The chemistry of the pnictides increases the cost of fiuc- 



tuations and renders the pnictides divalent, thereby elim- 
inating the quasiparticle multiplets. 

Our theory has several experimental consequences for 
the chalcogenides: there should be a quasiparticle dou- 
blet at positive energies 0.5 eV and 0.9 eV which can 
be probed by inverse photoemission techniques. Fur- 
thermore, the quasiparticle multiplets can be detected as 
side-peaks in the optical conductivity, again at 0.5 eV and 
0.9 eV. It would be interesting to study the temperature 
dependence of both the photoemission and optics near 
the coherence temperature (500 K) as we expect strong 
temperature dependence due to the many-body nature 
of the quasiparticle multiplets. Finally, the mechanism 
we outline is fairly general and applies to other corre- 
lated materials, provided the coherence temperature is 
large enough to be observed and less than the atomic 
multiplet splitting. 
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